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A derivation of the effective lagrangian for non-relativistic 
quantum chromodynamics and the heavy quarks effective field 
theory is given. Our calculation provides of a simple and 
systematic method of calculation of the full off shell effective 
lagrangian at tree level including all the 1/m corrections. 



Presently QCD is the best candidate for describing the 
strong interactions of its asymptotic freedom property at 
high energies. However, in the infrared domain new ad- 
ditional problems appear, namely, the effective coupling 
constant increase for low energies and as a consequence 
phenomena such as confinement, hadronization or chiral 
symmetry breaking must studied using non-perturbatives 
techniques. 

Non-relativistic quantum chromodynamics (NRQCD) 
and heavy quark effective field theory (HQET) are exam- 
ples of effective theories that describe approximately the 
low energy dynamics of QCD. Although both NRQCD 
|jj and HQET Q describe processes involving heavy 
quarks, they are physically different theories. Indeed, 
while NRQCD describes processes at low transferred mo- 
mentum, HQET only consider processes where kg >> 
k 2 /2m. From the physical point of view this difference 
is manifested in the structure of the quarks propagators, 
in NRQCD the propagator is 
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while in HQET is 
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However in spite of these physical differences, both the- 
ories are written in terms of an effective lagrangian com- 
ing from QCD in the non-relativistic limit. Normally this 
effective lagrangian is derived using symmetry arguments 
H and the coefficients of the operators involved in the 
expansion are obtained by using matching conditions [Q . 

In this letter we would like to propose a new ap- 
proach to the effective lagrangian calculation of NRQCD 



which has the advantage that it allows a derivation of 
all the relativistic corrections of the (tree-level) effective 
lagrangian, and in addition it is valid without using the 
equation of motion [||. 

In order to discuss our results let us start considering 
the QCD lagrangian 
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QCD 



Cq + C-L + C-H, 



(3) 



where Lq is the lagrangian for the gluons fields and Cl,h 
is the fermionic part for the light and heavy quarks re- 
spectively, i.e. 



£l,h = iPl,h[H? - m L .H\lpL.H, 



(4) 



where If) = ft - 

In order to define the heavy quark mass one must ne- 
glect the hard gluons contributions while the soft gluons 
one, by definition, are contained in the heavy quark mass. 
This assumption is valid when the heavy quark mass is 
much heavier than the scale of QCD and, under these 
conditions, the heavy quarks can be considered as non- 
relativistic particles. As a consequence of this, one focus 
the attention to the heavy modes sector of the partition 
function 



Z H [A] = / V$ H Vij) H e 



iS H 



(5) 



Heavy quarks interact with the light modes through 
the gluon field and this is weak if measured at the scale 
of the heavy fermion mass ( we will omit the subscript H 
from now on) . In this case the original bispinor ip can be 
written in terms of a slowly varying bispinor <fi as[] 



(6) 



where (f> , in the leading approximation, carries no infor- 
mation about the heavy quark mass. The only contribu- 
tion coming from the mass of the heavy quarks appears 
in the corrections in powers of 1/m. 

Using (jfy one find that the heavy quark lagrangian is 



£ = 4>{ilp - m(l - 7q))4>. 



(7) 



This can be written explicitly in terms of the large if and 
small x components of 4> as 
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1 Notice that this reparametrizations of the fields can be seen 
as changing the origin from where the energy E is measured, 
i.e., as defining Eq = E — m 
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£ = (pHD (p + x [iDo + 2m] \ + i& ■ D% 

+X ] ia-T) Vl (8) 

where we have used the standard Dirac's representation 
for gamma matrices. Thus, the partition function for the 
heavy quarks in terms of <p, and Xi is 



Z H [A] = J V^VipV^Vx e lSn . 



(9) 



The diagonalization of this lagrangian is straightforward. 
Indeed, if one performs the change of variables (with unit 
Jacobian) 



<P =<P, 
ip t = ip\ 

X = X + [iD + Zm}^ 1 ia ■ Dip, 
X ] = X f + ^ io ■ D [iD + 2m]" 1 



(10) 



in (|5[), the new lagrangian reads (omitting the primes) 

£ = ipi [iD Q + cr • D (iD Q + 2m)" 1 a ■ D] <p, 

+ X t [iD + 2m]x- (11) 

This lagrangian describes the (non-local) dynamics of 
relativistic heavy quarks in terms of two components 
spinors. One should note that ip and x appear decoupled 
and the integration in x contributes with a normaliza- 
tion factor (after choosing an appropriate gauge in this 
frame) . 

The non-relativistic limit of ([ll]) is straightforward be- 
cause one expands the operator (iD + 2m)" 1 in powers 
of l/m, i.e. 

{iDo + 2m)- 1 = -L(l - -i-Do + -^(Do? - ...). 
2m 2m 4m z 

(12) 

Then, the effective lagrangian for the heavy quarks com- 
ing from rtlll) becomes 



£ = £^ + 



(13) 



The first term in ( |l3| ) after to use OiOj — Sij + ieijk& k , 
reads 

=^t[ i£ , 0+ J_ D 2 + 5 . B ]^, (14) 
2m 2m 

which is just the lagrangian for a non-relativistic quark in 
a (chromo)magnetic field and it includes the Pauli term. 

Even if expanded in this way, the resulting lagrangian 
is not written in the standard form, i.e. £ — pq — Ti., from 
where one can infer directly the way the Hamiltonian 
looks like. To that end must use the equations of motion 
for both ip* and ip when calculating those contributions 
higher that £^ . 



The higher order corrections are more laborious to find 
but the calculation is straightforward. Thus, the order 
l/m 2 correction is 

= [[D, E] + ia ■ (D x E - E x D)] ip, (15) 

where the terms in the RHS are the well known non- 
abelian Darwin and spin-orbit ones, and in the conmuta- 
tor [D,E] an inner product is understood. 

The next higher order terms (l/m 3 ) can be computed 
following the same procedure, thus one gets that £^ is 



I g 2 a ■ (B x B - E x E)] ip, 



(16) 



where the second line is a genuine QCD contribution. 

The lagrangian dH) including 1/ m corrections was 
written using dimensional and plausibility arguments by 
Lepage et. al. in H . The above derivation provides us of 
a systematic calculation method where contact, spin and 
color terms simply does not exist in the non-relativistic 
limit. These last terms are discarded in the Lepage et. al. 
analysis || by energetic considerations and, as a conse- 
quence, the approach proposed here could be considered 
as a first principles calculation. 

The equation @ is the NRQCD la grangian in the 
rest frame, in an arbitrary frame (|l3|) becomes the HQET 
lagrangian. For details see B. 

We finalize this paper making some comments concern- 
ing to the change of variables (|Io[). The new variables 
diagonalize (||) producing a non-local term in (pd|). Al- 
though this non-local term retains all the information 
concerning to the relativistic corrections, one can see 
corrections such as Darwin, spin-orbit, etc. only order 
by order. Thus, we could conjecture that the operator 
(iDo + 2m) 1 contains all the information given by the 
Foldy-Wouthuysen transformation (|]. In fact, this trans- 
formation is applied to the equation of motion and then 
an extra condition is impossed in order to separate the 
small and large components of the spinor. In our case, 
things are done in reverse order, i.e., a transformation is 
obtained that makes the separation between small and 
large components and then the equations of motion are 
imposed. 

After this paper was completed , we were informed 
about reference H where a path integral derivation of 
the effective action for HQET was done . In a rest frame 
our formula ([ll]) corresponds to the Mannel et. al. ef- 
fective action (eq. (20) in this reference) . However we 
would like to emphasize the differences between both cal- 
culations; in our case instead of using projection opera- 
tors as in pL we propose the change of variables jl0| ) 
obtaining (O) straightforwardly. The explicit expansion 
of (iDo + 2m) 1 permit to derive eq. ( |T3| ) (i.e. (14)- 
(16)) instead of conjecturing it as it has been previously 
discussed. 
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